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ELASTOHYDRODYNAMIC FILM THICKNESS FORMULA BASED ON 
X-RAY MEASUREMENTS WITH A SYNTHETIC PARAFFINIC OIL 
by Stuart H. Loewenthal, Richard J. Parker, and Erwin V. Zaretsky 


Lewis Research Center and 
U. S. Army Air Mobility R&D Laboratory 

SUMMARY 

An empirical elastohydrodynamic (EHD) film thickness formula for heavily loaded 
contacts based upon X-ray film thickness measurements made with a synthetic paraffinic 
oil was developed. The empirical model was formulated from test data which covered a 
wide and practical range of operating conditions. Maximum Hertz stresses ranged from 
1. 04xl0 9 to 2. 42x1 0 9 N/m 2 (150 000 to 350 000 psi), disk temperatures from 339 to 
505 K (150° to 450° F), and mean surface speeds from 9. 4 to 37. 6 m/sec (370 to 1480 in. / 
sec). 

Predicted values of minimum film thickness from the deduced film thickness formu- 
la were compared to X-ray test data and to the numerical results from a theoretical iso- 
thermal EHD film thickness equation. In contrast to conventional theory, the present 
film thickness model did reflect the high sensitivity of film thickness to applied load ex- 
hibited by the test data. The empirical model showed improved agreement with the X-ray 
test data throughout the range of test conditions. 

It was judged that the empirical film thickness formula can be used with reasonable 
certainty for rolling-element bearing and gear systems composed of steel, employing the 
synthetic paraffinic oil studied herein and whose contact geometry approximates that upon 
which the model is based. 


INTRODUCTION 

The importance of maintaining a sufficient elastohydrodynamic (EHD) film thickness 
between dynamically contacting machine elements has in recent years been more fully 
appreciated. The prediction of elastohydrodynamic film thickness has been the focal 
point of many theoretical and experimental investigations (refs. 1 and 2). The ratio of 
the elastohydrodynamic minimum film thickness to the composite surface roughness has 



become an acceptable indicator of the effectiveness of the lubricant film within the rolling - 
element contact zone. It has been experimentally shown to influence the fatigue life of 
rolling- element bearings (refs. 3 and 4). Predetermination of this lubricant parameter 
involving an accurate prediction of minimum film thickness will be of great value to the . 
designer in obtaining realistic estimates of component fatigue life (ref. 5). 

The bulk of the experimental work conducted in elastohydrodynamic lubrication has 
been confined to conditions of moderate speeds; that is, up to 25. 4 m/sec (1000 in. /sec), 
and moderate loads; that is, maximum Hertz stresses to 1.24x10 N/m (180 000 psi) 
(refs. 6 to 9). The research of reference 10 has extended the elastohydrodynamic film 
thickness measurements to maximum Hertz stresses of 2.42xl0 9 N/m 2 (350 000 psi) 
which includes the design operating range of most machine components such as bearings 
and gears. In contrast to the results obtained by previous investigators which showed 
reasonably good correlation at moderate speeds and loads between elastohydrodynamic 
theory and film thickness measurement, the data of reference 10 showed a marked devia- 
tion between prediction and experiment. In particular, at high contact stresses; that is, 
maximum Hertz stresses greater than 1. 38xl0 9 N/m 2 (200 000 psi), the sensitivity of 
the film thickness to load as determined experimentally is far greater than that predicted 
by the classical elastohydrodynamic theory of references 11 and 12. 

Several attempts have been made to resolve the apparent discrepancy between theory 
and experiment. A critical examination of the X-ray technique itself was made (ref. 13) 
for possible load dependent experimental errors. However, no experimental factors 
were uncovered which could seriously alter the accuracy of the X-ray measurements. 

On the theoretical side, the influence of several possible rheological factors have been 
investigated, such as the effects of a non-Newtonian lubricant of the Ree-Eyring form 
(ref. 14), the effects of heating at the inlet of the contact region (ref. 15) and the effects 
of a i*educed lubricant viscosity -pres sure dependence using a composite exponential 
model (ref. 16) and using a power -law model (ref. 17). 

While each of the aforementioned modifications to elastohydrodynamic theory has 
succeeded somewhat in improving the agreement, between theory and experimental data 
within the heavy load regime, the resulting predicted values of film thickness differed 
little in magnitude from those computed using classical EHD theory. Furthermore, the 
modified theories do not sufficiently account for the high film thickness -load dependence 
to allow accurate predictions of film thickness under realistic operating conditions. 

It becomes the objective of the work reported herein to (1) develop an empirical 
elastohydrodynamic film thickness model based on an analysis of the experimental data 
of reference 10 and (2) compare the empirical relation derived to that of conventional 
elastohydrodynamic theory. 
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SYMBOLS 


a 

b 


CL 


E 1’ E 2 


E T 


H 


min 


h . 
mm 

k j 

n i 


n 


m 


Hz 


P H z 

Rjr R2 


R T 


T 

T 0 

u 


u 

u r u 2 

a 


minor semiaxis of Hertzian contact, m (in.) 
major semiaxis of Hertzian contact, m (in. ) 
coefficient defined indirectly by eq. (1) 
modulus of elasticity of elements 1 and 2, N/m 2 (psi) 




(psi) 


dimensionless minimum film thickness parameter, - mm 

R* 

film thickness in Cheng’s theory defined by eq. (8), m (in. 

minimum film thickness, m (in.) 

coefficient defined indirectly by eq. (6) 

exponent defined indirectly by eq. (1) 

exponent defined indirectly by eq. (3) 


p„ 

dimensionless stress parameter, ___ 

E T 

o 

maximum Hertz stress, N/m (psi) 

radius of elements 1 and 2 in rolling direction, m (in. ) 

(— + > m (in.) 

\ R i r 2 J 

exponent defined indirectly by eq. (6) 
local temperature, K (°F) 
disk temperature, K (°F) 

4_u 

dimensionless speed -viscosity parameter, 

E’R’ 

- (uj + u 2 ), m/sec (in. /sec) 

2 


surface velocities of elements 1 and 2, m/sec (in. /sec) 

2 ? 

pres sure -viscosity coefficient, m /N (in. /lbf) 



2 2 

contact inlet absolute viscosity, N- sec/m , (lb-sec/in. ) 

Poisson T s ratio of elements 1 and 2 

EHD FILM THICKNESS EXPERIMENTAL DATA 

The film thickness model which is presented herein was developed from film thick- 
ness data (ref. 10) acquired in an X-ray rolling disk machine (fig. 1) with a synthetic 
paraffinic oil. This particular oil is of interest since it is one of a class of several ad- 

r Container-heater unit 


\ 

\ 



Figure 1. - X-ray rolling-contact disk machine (ref. 10). 


vanced lubricants singled out for high temperature jet engine application. Its performance 
as a bearing lubricant has been extensively evaluated (e. g. , see refs. 18 and 19). 

The method of measuring film thickness with the X-ray technique comprises pro- 
jecting X-rays between the surfaces of the two contacting disks and detecting the rate of 
X-ray transmission through the contact. Since the greatest constriction occurs at the 
trailing edge of the contact, the X-ray count thus becomes a measure of the lubricant's 
minimum film thickness. 
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The range of test conditions include disk temperatures from 339 to 505 K (150° to 

450° F), surface speeds from 9.4 to 37.6 m/sec (370 to 1480 in. /sec) corresponding to 

disk rolling speeds from 5000 to 20 000 rpm, and maximum Hertz stresses from 1.04x10® 

to 2.42x10® N./m^ (150 000 to 350 000 psi). Two crowned-cone AISI M-50 steel disks 

each with a radius of 1. 83 centimeters (0. 72 in. ) and a surface finish of 2. 5x10"^ to 
6 

5.0x10'° centimeter (1 to 2 p in.) rms were used as the test specimens (fig. 2). (Both 
crowned disks and crowned-cone disks (with a cone angle of 10°) were tested and no sig- 
nificant differences were found between the two sets of data. ) All the data reported here- 
in were generated with the crowned-cone test disks. The properties of the synthetic par- 
affinic oil which was used in the course of this study are listed in table I. 



Figure 2. - Contacting disk geometry (ref. 10). (All linear dimensions 
in cm fin. ).) 


TABLE I. - PROPERTIES OF THE SYNTHETIC PARAFFINIC OIL 


Kinematic viscosity, cS (or 10”^ m^/sec), at - 

233 K (-40° F) >100 000 

311 K (100° F) 443 

372 K (210° F) 40 

478 K (400° F) 5.8 

589 K (600° F) a 2.3 

Flash point, K (°F) 541 (515) 

Fire point, K (°F) 589 (600) 

Autoignition temperature, K (°F) 703 (805) 

Volatility (6.5 hr at 533 K (500° F)), wt. % 14.2 

Specific heat at 533 K (500° F), J/(kg)(K) (Btu/(lb)(°F)) 2910 (0. 695) 

Thermal conductivity at 533 K (500° F), J/(m)(sec)(K) (Btu/(hr)(ft)(°F)) 0.12 (0.070) 

Specific gravity at 533 K (500° F) 0.71 

Pressure-viscosity coefficient at 298 K (77° F), m^/N (psi' 1 ) b 1.3xl0 - ^ (0. 92x10'“*) 


a Extrapolated. 
^From ref. 25. 
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EHD FILM THICKNESS MODEL 


A summary of the X-ray test results showing measured minimum film thickness 
h min P lotted against maximum Hertz stress p Hz for various mean surface speeds u 
and disk temperatures Tq is presented in figure 3. Within a heavily loaded contact, the 
film thickness dependence upon load as evidenced by the X-ray data is far greater than 
any current EHD theory predicts. This is illustrated by the dashed line appearing in 
figure 3(a) whose slope represents the typical stress exponent from a commonly used 



I I ^ 1 1 1 1 5 

1.5 2 3 4X10 5 1.5 2 3 4x10° 

Maximum Hertz stress, p Hz , psi 


Disk speed, 
rpm 

O 5 000 
□ 10 000 
A 15 000 
O 20 000 



l I l I I , 

1.5 2 3 4xl0 5 

Maximum Hertz stress, p Hz , psi 


(cl Disk temperature, 422 K 
GOO 0 F). 


(d) Disk temperature, 478 K 
(400° FI. 


Ie) Disk temperature, 505 K 
(450° F). 
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Figure 3. - X-ray test results showing effect of maximum Hertz stress on measured minimum film thickness with a synthetic paraffinic 
oil using crowned-cone disks tfram ref. 7). 



film thickness formula (ref. 16). The pronounced effect of increasing load on h . is 
clearly apparent from these graphs. 

By introducing the customary nondimen sional elastohydrodynamic groupings, the 
measured X-ray data can be represented over the entire range of test conditions in a 
single plot. This is accomplished in figure 4, which is a log-log plot of the dimension- 
less film thickness parameter H m ^ = ^min^®* against the dimensionless speed- 

viscosity parameter U = /i u/E T R’ for various maximum Hertz stresses. The straight 


Maximum Hertz stress, 
N/m 2 <psi> 

o 1.04xI0 9 <150 000) 

--D-- 1.37 ( 200 000) 

0 1.72 (250 000) 

--a- — 2.07 ( 300 000) 



Dimensionless speed-viscosity parameter, U 

E'R 


Figure 4. - Effect of maximum Hertz stress on sensitivity of minimum film thickness to changes in speed and 
viscosity for synthetic paraffinic oil. 
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lines appearing on this figure have been least-square fitted through the measured data for 
each test stress level and can be represented by the following simple power relation: 



or in dimensionless groupings 


H min 



( 1 ) 


( 2 ) 


where parameters C. and n., listed in table H, are some function of maximum Hertzian 
contact pressure, say f(p Hz ). 

It is apparent from figure 5 that the measured minimum film thickness shows an 
enhanced sensitivity to surface speed and lubricant viscosity with increasing stress level. 
The variation of exponent m with operating conditions is not accounted for by conven- 
tional elastohydrodynamic theory. However, Westlake and Cameron (ref. 20), in studying 
some 35 different fluids, report a similar variation of the exponent for the nondimension - 
al speed-viscosity parameter with contact load for their empirical midfilm thickness 
equation. These data were generated using an optical interferometric method for deter- 
mining the amount of film generated between a rolling steel ball and a flat glass plate. 

It should be mentioned that no such variation of exponent with stress level was observed 
for their minimum film thickness data. Furthermore, the data of reference 20 was ob- 
tained from elastohydrodynamic contacts in which the lubricant was subject to shear rate 
levels and contact pressure levels which are considerably lower than those generated by 
the rolling disk contact in the X-ray experiments. 

For purposes of developing an approximate film thickness expression for synthetic 
paraffinic oils, the complication of a pressure dependent speed -viscosity parameter ex- 
ponent can be avoided without introducing serious inaccuracies by selecting some mean 
value of n., that is, n m , so that equation (2) can be written as 


H min 


n 

u ■ 


f <0Hz> 


(3) 


where the function f(Pjj z ) has been introduced to mathematically model the effect of con- 
tact stress on minimum film thickness and where it has been determined that n m = 0. 66. 
By comparing equations (2) and (3), it is clear that 


_(n.-0. 66) 

f(p Hz ) « c 4 u 


(4) 


8 



TABLE II. - EMPIRICAL MINIMUM FILM THICKNESS COEFFICIENTS FOR EQUATION (1) 


Dimensionless 

stress 

parameter, 

p _ P Hz 

P Hz -^T 

Maximum Hertz stress 

Coefficient 

C i 

Exponent 

n. 

l 

N/m 2 

psi 

3.09xl0" 3 

1. 04xl0 9 

150 000 

28. 1 

0. 584 

4. 12 

1. 38 

200 000 

44.6 

. 609 

5. 15 

1. 72 

250 000 

93.5 

.649 

6. 17 

2.07 

300 000 

231. 

. 702 

7.20 

2.42 

350 000 

654. 

.764 



Figure 5. - Sensitivity of minimum film thickness 
to maximum Hertz stress as a function of surface 
speed and lubricant viscosity for synthetic paraf- 
finic oil. 




or more concisely 


f kW = 


H min 

u °- 66 


(5) 


It is evident from equation (4) that f(p Rz )_ is unavoidably an explicit function of U. 
The extent of the dependence of f(p Rz ) with U over its experimental range is illustrated 
in figure 5. In this plot the broken lines represent the variation of the log of f(p Rz ) with 
the log of the dimensionless stress parameter P Rz at the minimum and maximum exper- 
imental values of U. It is apparent from this figure that the sensitivity of minimum film 
thickness to maximum Hertz stress will vary slightly in accordance with the operating 
speed and temperature condition. At the risk of introducing a small degree of error, the 
extent of which will become evident later, and in order to conform more closely to con- 
ventional elastohydrodynamic theory, the variation of film thickness to contact stress will 
be modeled for a constant mean value of the experimental dimensionless speed-viscosity 
parameter. This variation is represented by the solid curve in figure 5. There are^ob- 

viously several ways of mathematically representing this variation of f(Pjj z ) w ith ^Hz‘ 

The one which is presented here is partially in keeping wiih conventional eiasiuhycuoay- 
namic theory, that is, let 


f(P HZ > = Vhz S] (6) 

where coefficients K. and S. would normally be constants independent of contact load 
under conventional theory but^ere by necessity K- and S- will be deemed as quasi- 
constant, that is, constant over specific intervals of operating contact stress as shown 
in table HI. Thus when equations (3) and (6) are combined, they form a complete cor re - 


TABLE TIL - MINIMUM FILM THICKNESS CORRELATION COEFFICIENTS FOR EQUATION 7 


Range of operating contact stress 

Coefficient 

K i 

Exponent 

s i 

Dimensionless stress 

P'U'r 

parameter, P Hz = 

E' 

Maximum Hertz stress 

. 9 
N/m 

psi 

3.09xl0' 3 to 4. 12xl0' 3 

1 . 04x1 0 9 to 1. 38xl0 9 

150 000 to 200 000 

31. 1 

- . 38 

4. I3xl0" 3 to 5. 15x10" 3 

1. 39x1 0 9 to 1. 72xl0 9 

201 000 to 250 000 

2.81 

- . 72 

5. 16xl0’ 3 to 6. I7xl0" 3 

1 . 73xl0 9 to 2.07xl0 9 

251 000 to 300 000 

.0864 

-1. 38 

6. 18x10’ 3 to 7. 20x10” 3 

2- 08xI0 9 to 2. 42xl0 9 

301 000 to 350 000 

.00127 

-2.21 


10 




lation for determining minimum film thickness over a wide range of temperature and 
speed conditions for heavy contact loads, that is, above a Pg z of 1.04xl0 9 N/m (150 000 
psi). This correlation can be written: 


H . = K.U°' 66 P H7 S;j 
min j Hz 


( 7 ) 


where the four discrete values of coefficients and S- over the maximum Hertz stress 
range from 1. 04xl0 9 to 2. 42xl0 9 N/m 2 (150 000 to 350 000 psi) appear in table HI. 


DISCUSSION OF RESULTS 


Comparison with Conventional EHD Theory 


A comparison of the film thickness correlation from equation (7) with the experi- 
mental X-ray data (ref. 10) is shown in tabular form in table IV and in graphical form 
with the isothermal film thickness equation of Cheng (ref. 21) in figures 6, 7, and 8. 
Cheng’s film thickness formula is considered typical of those from conventional elasto- 
hydrodynamic theory for line contact. It can be written 


h 

-- = 1.47 
R’ 



( 8 ) 


In computing Cheng’s film thickness, a correction factor of 0. 8 has been applied to 
adjust Cheng’s center film thickness to minimum film thickness. Values of the pressure- 
viscosity coefficient a, for the synthetic paraffinic oil utilized in this computation were 
extracted from reference 22 and appear in table V as a function of temperature. In this 
report (ref. 22) several parameters representing the pressure -viscosity characteristics 
for this oil were presented. The a parameter deemed most appropriate for film thick- 
ness calculations to be made here was reportedly obtained by applying a numerical inte- 
gration technique to the experimental data. The necessity of extrapolating the published 
data for a over a relatively large temperature interval (83 K (150° F)), places the value 
estimated in table V for a at 505 K (450° F) in some question. 

As anticipated, it is apparent from the plots presented in figures 6, 7, and 8 that 
there is improvement in predicting actual minimum film thicknesses with the present 
empirical model (eq. (7)) for the system under study. This improvement is largely 
attributed to the ability of the deduced film thickness relation to reflect the high sensi- 
tivity of minimum film thickness to contact stress under heavy applied loads. It is evi- 
dent that deviations of up to several hundred percent exist between measured and pre- 
dicted values of minimum film thickness using commonly accepted EHD formulas. 
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Dimensionless film thickness parameter, 


Cheng's formula (eq. (8)1 

Correlation leq. (71) 

X-ray test data (from fig, 3) 


500xl0' 6 



(a) Disk speed, 51)00 rpm. (Q I Disk speed, 10 000 rpm. 

9D0x 10" 6 

Km P- f— 



Dimensionless stress parameter, P Hz = 

(c) Disk speed, 15 000 rpm. (d) Disk speed, 20 000 rpm. 

Figure 6. - Comparison of predicted minimum film thickness 
with X-ray test data for synthetic paraffinic oil at disk temper- 
ature i r 0 - 339 K (150° F). 


Cheng's formula (eq. (8)) 

Correlation (eq. (7)) 

X-ray test data (from fig. 3) 


Afixin' 6 



Dimensionless stress parameter, P^ z = 

(c) Disk speed, 15 000 rpm. (dl Disk speed, 20 000 rpm. 

Figure 7. - Comparison of predicted minimum film thickness with 
X-ray test data for synthetic paraffinic oil at disk temperature 
Tq = 422 K (300° FI. 





TABLE IV. - COMPARISON OF EMPIRICAL FILM THICKNESS MODEL WITH X-RAY DATA 


FOR A SYNTHETIC PARAFFINIC OIL (REF. 7) 


(a) SI units 


Disk 

tem- 

pera- 

ture, 

K 

Disk 

speed, 

rpm 

Maxi- 

mum 

Hertz 

stress, 

N/m 2 

Minimum film 
thickness, 

h min' cm 

Disk 

tem- 

pera- 

ture, 

K 

Disk 

speed, 

rpm 

Maxi- 

mum 

Hertz 

stress, 

N/m 2 

Minimum film 

thickness 

h . , p cm 
min 

Disk 

tem- 

pera- 

ture, 

K 

Disk 

speed, 

rpm 

Maxi- 

mum 

Hertz 

stress, 

N/m 2 

Minimum film 
thickness, 

h min> » cm 

X-ray 

Corre- 

lation 

X-Ray 

Corre- 

lation 

X-Ray 

Corre- 

lation 

330 

5 000 

1.04xl0 9 

89 

74 

422 



20 

15 

505 


1.04x10° 

5 

8 



1. 38 

86 

69 



1. 38 

18 

15 



1. 38 

5 

8 



1. 72 

71 

56 



1. 72 

10 

13 



1.72 

5 

5 



mm 

El 

46 



2,07 

5 

10 




2 

5 



BE 

mb 




2.42 

5 

8 



2.42 

2 

2 



1.04xl0 9 




10 000 

1.04x10° 

33 

25 




12 

13 



1.38 

99 




1. 38 

30 

23 



1. 38 

10 

10 



1. 72 

86 

BK 



1. 72 

20 

20 



1.72 

a 

10 



2.07 

WM 

■sit 



2.07 

12 

15 



2.07 

5 

8 



2.42 

■9 

KB 



2.42 

10 

10 



2.42 

2 

5 



1.04x10° 

119 

150 


15 000 

1.04x10° 

41 

33 




15 

Dl 



1. 38 

109 

140 



1. 38 

35 

30 



1. 38 

15 

Hi 



1. 72 

94 

119 



1.72 

25 

25 



l. 72 

10 

H9 



2.07 

81 

91 



2.07 

15 

20 



2.07 

8 

Bf 



IESH 

53 

66 



2.42 

15 

15 



2.42 

5 

Rfl 



1.04x10° 

127 

EX 


20 000 

1.04x10° 

43 

41 



1.04x10° 

18 




1.38 

117 




1.38 

41 

38 



1. 38 

18 

18 



1. 72 

102 




1 72 

30 

30 



1.72 

15 

15 



2.07 

84 




2.07 

20 

25 



2.07 

10 

10 



2.42 

58 




2.42 

20 

18 



2.42 

8 

8 

366 

5 000 

1.04x10° 

43 


478 

5 000 

1.04x10° 

12 

10 








1. 38 

41 




1.38 

10 

8 








1. 72 

38 




1.72 

5 

8 








2. 07 

25 




2.07 

2 

5 








2.42 

20 




2.42 

2 

5 







10 000 

1. 04x10° 

61 



10 000 

1.04x10° 

23 

15 








1.38 

61 




1. 38 

23 

13 








1. 72 

58 





8 

13 








2.07 

38 

38 




10 

10 








2.42 

30 

25 



1 

8 

8 







15 000 

1.04x10° 

66 

mi 


15 000 


30 

20 








1. 38 

66 




1. 38 

28 

18 








1. 72 

64 




1. 72 

18 

15 








2.07 

43 

PH 



2.07 

15 

13 








2.42 

38 




2. 42 

10 

ft 







20 000 

1.04x10° 

71 



20 000 

1.04x10° 

36 

23 








1. 38 

71 




1. 38 

30 

20 








1. 72 

66 

K9 



1. 72 

20 

18 








2.07 

48 

si 



2.07 

18 

15 









2. 42 

41 

KB 



2.42 

15 

10 
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TABLE IV. - Concluded. COMPARISON OF EMPIRICAL FILM THICKNESS MODEL WITH X-RAY DATA 


FOR A SYNTHETIC PARAFFINIC OIL (REF. 7) 


(b) U.S. customary units 


Disk 

tem- 

pera- 

ture, 

°F 

Disk 

speed, 

rpm 

Maxi- 

mum 

Hertz 

stress, 

psi 

Minimum film 
thickness, 

h min> ^ in ’ 

Disk 

tem- 

pera- 

ture, 

°F 

Disk 

speed, 

rpm 

Maxi- 

mum 

Hertz 

stress, 

psi 

Minimum film 

thickness 

h . , p cm 
min 

Disk 

tem- 

pera- 

Disk 

speed 

rpm 

Maxi- 

mum 

Hertz 

stress, 

psi 

Minimum film 
thickness, 
h min' ycm 

X-ray 

Corre- 

lation 

X-Ray 

Corre- 

lation 

°F 

X-Ray 

Corre- 

lation 

150 


PfiPfTiy 

35 


300 

5 000 

1.5xl0 5 

8 

6 

450 



B 

IB 




34 




2.0 

7 

6 


mm 


wm 

Hi 



2.5 

28 

23 



2.5 

4 

5 




MS 





23 

18 



3.0 

2 

4 




H 




3.5 

16 

13 



3.5 

2 

3 


i * 

Ey9 

■ 

'SH 


10 OOO 


42 

45 


10 000 

1.5xl0 5 

13 

10 



1.5x10 5 


5 




39 

42 



2.0 

12 

9 





4 



2.5 

34 

36 



2.5 

8 

S 



2. 5 


4 




28 

28 



3.0 

5 

6 



3.0 


3 


m m 1 

3.5 


20 



3.5 

4 

4 



3.5 


2 



1. 5xl0 S 


59 


15 000 

I.5xlO S 

16 

13 



1. 5xl0 5 

6 

6 



2.0 


55 



2.0 

14 

12 



2.0 

6 

6 



2.5 


47 



2-5 

10 

10 



2.5 

4 

5 



3.0 


36 



3.0 

6 

8 



3.0 

3 

4 



3.5 


26 



3.5 

6 

6 



3. 5 

2 

3 


20 000 


50 

71 


20 000 

1. 5xl0 5 

17 

16 


20 000 

1. 5xl0 5 


7 



2.0 

46 

66 



2.0 

16 

15 



2.0 





2.5 

40 

56 



2.5 

12 

12 



2.5 





3.0 

33 

44 



3.0 

8 

10 



3.0 





3.5 

23 

31 



3.5 

8 

7 



3.5 



200 

5 000 

1.5xl0 5 

17 

15 

400 

5 000 

1. 5x10° 

5 

4 

H 







2.0 

16 

14 



2.0 

4 

3 








2.5 

15 

12 



2.5 

2 

3 








3.0 

10 

9 



3.0 

mm 

2 








3.5 

8 

7 



3.5 

n 

2 








1.5x10® 

24 

24 


10 000 

1.5xl0 5 

H 

G 








2.0 

24 

22 



2.0 

■9 

5 








2.5 

22 

19 



2.5 

Kg 

5 








3.0 

15 

15 



3.0 

H 

mm 








3. 5 

12 

10 



3.5 

3 

Isi 







15 000 


26 

31 


15 000 

1.5xl0 5 

12 

m 








2.0 

26 

29 



2.0 

E9 

mm 








2.5 

25 

25 



2. 5 

19 









3.0 

17 

19 



3.0 

6 

wm 








3.5 

15 

14 



3.5 

4 

wm 







20 000 

1.5xl0 5 

28 

38 


20 000 

1.5xl0 5 



■ 







2.0 

28 

35 



2.0 


8 








2.5 

26 

liiSS 



2.5 


7 








3.0 





3.0 


6 








3.5 

16 

' 1 



3.5 

6 

4 
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TABLE V. - ESTIMATED VALUES OF PRESSURE -VISCOSITY COEFFICIENT o 


FOR A SYNTHETIC PARAFFINIC OIL (FROM REF. 22) 


Disk 

temperature 

Pressure -viscosity 
coefficient, 
a 

K 

°F 

m^/N 

psi -1 

339 

150 

17. lxlO -9 

11.8xl0' 5 

422 

300 

10.8 

7. 5 

505 

450 

a 4.9 

a 3. 4 


a Extrapolated. 


Application of the empirical film thickness expression to systems where different lubri- 
cation conditions prevail will be discussed next. 


Effects of Contact Geometry, Material, and Lubricant 

The present correlation is based exclusively upon measurements made with a single 
test lubricant and disk geometry and material combination. The disk geometry was 
chosen to simulate the ball-inner race contact of a 120-millimeter bore angular contact 
ball bearing (ref. 23). The contact between the test disks approaches the condition of 
line contact with an ellipticity ratio b/a of 5. 9 where b and a are, respectively, the 
major and minor semi-axes of the contact ellipse. The equivalent radius of curvature 
in the direction of rolling R* for the test disks was 0. 915 centimeter (0. 36 in). 

It is difficult to extend the results presented herein to different contact geometries, 
materials, and lubricants without further experimental verification. Indeed, each of 
these factors are known to influence film thickness in varying degrees. For example, 
Cheng (ref. 21) has theoretically shown that the proportions of the contact ellipse with 
b/a varying from 1 to 5 have a relatively mild effect on film thickness. The dependence 
of film thickness upon mean surface speed u, absolute viscosity and contact stress 
pjj z changes little as the shape of the contact ellipse varies from point to line contact. 
Similarly, Archard and Cowking in reference 24 have shown that there is great similar- 
ity between the EHD lubrication of point and line contacts. In addition, the typical con- 
tacts between the races and the balls or rollers in rolling -element bearings and the con- 
tacts between gear teeth approximate the line contact case. In view of the aforemen- 
tioned, the minimum film thickness correlation (eq. (7)) is suitable for most practical 
applications without further modifications for small differences in contact geometry. 
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With regard to the size of the contacting elements, elastohydrodynamic theory indi- 
cates that film thickness is moderately dependent upon the contacting element’s equiva- 
lent radius of curvature in the rolling direction R'. In utilizing nondimen sional group- 
ings, the film thickness correlation implies that film thickness is a function of R T to the 
0. 34 power at a given contact stress level. This is reasonably in accord with what con- 
ventional EHD theory forecasts; for example, see Cheng’s formula (eq. (8)) where 

0 9fi 

h Q ^R for p Hz = constant. The sensitivity of minimum film thickness to the con- 
tracting element's R’ in the heavy load regime remains to be established experimentally. 
Until such time, minimum film thicknesses predicted by the deduced relation in equa- 
tion (7) will be most meaningful for those systems in which the R' of the contact approx- 
imates the R' of the test disks. \ 

The effect of the material parameter E’, which is termed as the equivalent modu- 
lus of elasticity, is reflected in the film thickness expression through the dimensionless 
pressure parameter P Hz and the dimensionless speed-viscosity parameter U. But due 
to the varying nature of the exponent of P^ z , the effect on calculated film thickness of 
materials having elastic properties different from steel would be unjustifiably great at 
high contact loads. Hence, the film thickness correlation presented here should only be 
used for those rolling-element bearing and gear systems composed of steel 
(E* = 33. 6xl0 10 N/m 2 (48.6x10® psi)). 

It is expected that the most significant effect on equation (7) will be due to lubricant 
properties. The empirical coefficients utilized in the correlation, namely, n m , Ky 
and S-, have been derived from the experience gained with a single test fluid. Conse- 
quently, no meaningful generalizations can be made with regard to the effects of different 
lubricant types or formulations on the values of these coefficients without benefit of addi- 
tional experimental data. Until such time, the application of the film thickness correla- 
tion must be restricted to those systems employing the synthetic paraffinic oil upon which 
the deduced relationship is based. 

It may be apparent that the pressure-viscosity coefficient a which is customarily 
used to characterize the film forming capabilities of a lubricant apart from the effects of 
the lubricant’s absolute viscosity, is conspicuously absent from the present correlation. 

In the present model, the role formerly played by the pressure-viscosity coefficient has 
been fulfilled in part by the lubrication coefficient K-. That is, the relative value of a 
lubricant's K- and the value of its absolute viscosity fi Q at given operating condition is 
indicative of the amount of film thickness this particular lubricant can be expected to 
generate. An important distinction between a and K., is that coefficient a is solely 
temperature dependent while coefficient is solely contact stress dependent. Careful 
examination of the X-ray data revealed that the effect of temperature on film thickness 
could be adequately reflected by the /i Q variable in the dimensionless speed -viscosity 
parameter u, and thus the added complication of an additional temperature dependent 
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variable could be avoided. Secondly, there has been some difficulty in properly charac- 
terizing the pres sure -viscosity characteristics of fluids for purposes of calculating film 
thickness. This is due mainly to difficulties one encounters when attempting to develop 
a test apparatus which can accurately reproduce the very severe operating pressure and 
shear stress levels which a lubricant experiences within an elastohydrodynamic contact 
(ref. 20). Consequently, relevant test data of this nature has been generally limited. In 
view of the aforesaid, it was judged advantageous to dispense with ex for the present 
film thickness model. 

In conclusion, the film thickness relation developed herein represents an initial at- 

9 2 

tempt at empirically modeling the effects of high contact stress (above 1.04x10 N/m 
(150 000 psi) on minimum film thickness in an elastohydrodynamic contact. It is hoped 
that this formula, which does reflect the high load dependence exhibited by experimental 
data, will aid the designer in obtaining a better appraisal of the extent of lubricant film 
separating his contacting machine elements. As more experimental data becomes avail- 
able, the film thickness model will be extended to other lubricants and different operating 
parameters. 


SUMMARY OF RESULTS 

An empirical elastohydrodynamic (EHD) film thickness formula for heavily loaded 
contacts based upon X-ray film thickness measurements made with a synthetic paraffinic 
oil was developed. The minimum film thickness formula was developed from test data 
which covered a wide and practical range of operating conditions. Maximum Hertz 
stresses ranged from 1. 04x10^ to 2. 42x10^ N/m^ (150 000 to 350 000 psi), disk tempera- 
tures from 339 to 505 K (150° to 450° F), and mean surface speeds from 9. 4 to 37. 6 
m/sec (370 to 1480 in. /sec). Predicted values of minimum film thickness were com- 
pared to the X-ray test data for the same conditions and to the results from a theoretical 
isothermal EHD analysis. The effects of contact geometry, material, and lubricant 
factors upon predicted film thickness were considered. The following results were ob- 
tained: 

1. In contrast to present elastohydrodynamic theory, the empirical film thickness 
model does reflect the high sensitivity of minimum film thickness to contact stress ex- 
hibited by the test data under heavy applied loads. Deviations of up to several hundred 
percent were observed between measured and predicted values of minimum film thick- 
ness and those forcasted by a typical film thickness formula at high contact stresses. 

2. The measured minimum film thickness data were observed to display an enhanced 
sensitivity to mean surface speed and lubricant absolute viscosity with increasing con- 
tact stress. 
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3. Finally, it was judged that the empirical film thickness formula can be used with 
reasonable certainty for rolling- element bearing and gear systems composed of steel, 
employing the synthetic paraffinic oil specified herein, and whose contact geometry 
approximates that upon which the model is based. 

Lewis Research Center, 

National Aeronautics and Space Administration, 
and 

U. S. Army Air Mobility R&D Laboratory, 

Cleveland, Ohio, August 7, 1973, 

501-24. 
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